Electronic properties of alkali-metal loaded zeolites 
— a "supercrystal" Mott insulator 
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First-principles band calculations are performed for the first time for an open-structured zeolite 
(LTA) with guest atoms (potassium) introduced in their cages. A surprisingly simple band structure 
emerges, which indicates that this system may be regarded as a "supercrystal" , where each cluster of 
guest atoms with diameter ~10A acts as a "superatom" with well-defined s- andp-like orbitals, which 
in turn form the bands around the Fermi energy. The calculated Coulomb and exchange energies for 
these states turn out to be in the strongly-correlated regime. With the dynamical mean-field theory 
we show the system should be on the Mott-insulator side, and, on a magnetic phase diagram for 
degenerate-orbital systems, around the ferromagnetic regime, in accord with experimental results. 
We envisage this class of systems can provide a new avenue for materials design. 

PACS numbers: 82.75.Vx, 71.30.-fh, 36.20.Kd 
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Materials design is usually done by varying chemical 
elements and/or atomic crystal structures. Here we con- 
sider doped zeolites, which are unique in that (i) they are 
host-guest systems, where the host itself appears in a rich 
variety (crystal structure of cages, size of the inter-cage 
window, etc), on top of rich possibilities for the species 
and the number of guest atoms [1]. (ii) the size of the 
cage is ^ nanometer. We can note the distinction of 
the zeolites from the solid fuUerene: The latter are also 
nanostructured supercrystals, but the relevant orbitals, 
even when doped with alkali metals, are LUMO/HOMO 
of the cage (buckeyball) rather than those of dopants [2] . 
While (i) implies higher degrees of freedom for the design, 
(ii) implies that the Coulomb interaction energy for the 
cluster doped in the cage, which only decreases inversely 
with the size, can remain large (~ eV for nanometer di- 
mensions). Here we envisage that these can provide a 
fascinating playing ground for systematically control the 
electron correlation as a new avenue for materials design. 

Indeed, experimental results by Nozue et al showed 
some zeolites (Fig.l) loaded with clusters (~ five atoms 
per cage) of potassium is ferromagnetic for T < 8K [3,4], 
which is unique in that magnetism occurs even though 
all the ingredients are non-magnetic elements [5] . At the 
same time, an infrared analysis for the K-loaded zeolite 
shows clear absorption edges, indicative of insulating be- 
havior [6]. 

Here we address both phenomena, and show they can 
be explained in terms of strong electron correlations in 
a surprisingly simple electronic structure. We start by 
showing that the relevant states of K-doped zeolite can be 
viewed from a "supercrystal" point of view. Namely, the 



single-particle electronic structure as computed within 
the local-density approximation (LDA) for a typical ze- 
olite having a cubic array of cages is very accurately de- 
scribed as a tight-binding band of "superatoms" (i.e., the 
states in the clusters of dopants) with simple s and p or- 
bitals. This explains why the optical spectrum can be 
interpreted as transitions between s and p (or p and d) 
orbitals in a simplified well [7] . For the LDA calculation 
we adopt the all-electron full-potential linear muffin-tin 
orbitals (FP-LMTO) [8]. In fact ours is the first reliable 
band structure calculation for zeolites [9]. 

Next we estimate the Coulomb and exchange interac- 
tions for these orbitals by plugging in the LDA wave func- 
tions, and show that the interactions are large compared 
with the single-particle bandwidth, indicating strong 
electron correlations. Last, we study the effect of cor- 
relations on the spectral function in the spirit of the 
LDA-|-dynamical mean- field theory [10], and show the 
system resides well on the Mott-insulator side. This re- 
solves a puzzle that the band calculation predicts that 
the system is a metal, while experiments indicate an in- 
sulator; moreover we show the calculated exchange en- 
ergy (Hund's coupling in superatoms) indicates ferromag- 
netism in the multi-orbital magnetic phase diagram, in 
accord with experiments. 

We start with the undoped zeolite with a typical struc- 
ture having a simple-cubic array of cages (Fig.l), where 
each cage (called a) is an Archimedes polyhedron. The 
region surrounded by eight a cages forms another cage 
called (3. Aluminosilicate with this structure (zeolite 
A) has Si and Al atoms situated at the vertices of the 
cages and connected by O atoms with the inner di- 
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ameter of a cage 10 A. The material used in most 
experiments [3,4] is K(potassium)-form zeolite A (ab- 
breviated as LTA hereafter), whose chemical formula, 
Ki2Ali2Sii2048 (with 84 atoms per imit cell), already 
contains some K atoms. Eleven out of the twelve K atoms 
are located on the centers of faces of the a cage (orange 
atoms in Fig.l), while the remaining K (red), assumed 
to be at the cage center, changes its position as K's are 
added to change the doping as we shall describe. For the 
atomic positions of Al, Si, O within the unit cell, we have 
adopted the accurate crystal structure obtained from a 
recent neutron powder diffraction study by Ikeda et al 
[11,12]. 

The result for the band structure for the undoped LTA 
is displayed in Fig.2, which shows that the undoped ze- 
olite is an insulator with a gap ~ 4 eV, where the wave 
functions indicate that, while the valence-top state sits 
on the framework, the conduction-band states (labeled 
as B-F) primarily reside within the a cage (for B,D,E,F) 
or within /3(C). 

The dispersion of the lowest conduction band (B) can 
be fitted, as shown in the right panel of Fig.2, by the 
dispersion of a tight-binding band on the simple cubic 
lattice, e(k) = txCos{kx) +tyCOti{ky) + tzCos{kz), with 
itx,ty,t,_) ~ (-20,-10,-5) (meV). The fitting is excel- 
lent, and, with all of tx, ty, < and the wave function 
being nodeless (Fig. 2(B), with its tetrahedral shape due 
to the configuration of surrounding K atoms), we may 
interpret the band as an "s band" of the supercrystal. 
Similarly, three (D, E, F) out of the four next-lowest con- 
duction bands, with large amplitudes within the a cage, 
may be interpreted as p bands. 

Now we come to the doped case. When one K atom 
is doped per unit cell (denote as KiLTA), the Fermi 
level shifts to the middle of the conduction (s) band 
(not shown), while the p band remains empty. Tight- 
binding fits to these bands arc again excellent, with 
{tx,ty,tz) — (—10,— 10,— 5) in mcV for the s band and 
~ (25,-10,-1) for the p. The good fit to the simple 
tight-binding model is highly nontrivial, since the frame- 
work, (K+)i2(AP+)i2(Si''+)i2 02")48 with nominal va- 
lence indicated, is an ionic compound. Chemically, the 
cage has a low electron affinity so that the electrons stay 
well away from the wall. This is reflected in the local- 
ization of the wave function, and provides an intuitive 
reason why a simple well is effectively realized. 

If we further dope the system to have K3LTA (which 
is experimentally in the magnetic regime), we have now 
four K atoms (the doped three on top of a red ball in 
Fig. 1(b)), which form a cluster in the a cage. The pre- 
cise atomic configuration has not been experimentally 
determined, so we have focused here on the configura- 
tion where the four atoms (at (0.5 ± a;i,0.5 ± a;i,0.5)) 
form a square, with ,ti ^ 0.25 minimizing the total en- 
ergy, although we have calculated for other possibilities 
such as a tetrahedral configuration. 



If we look at the band structure and the wave func- 
tions for the K3LTA in Fig. 3, we can see that three 
bands around Ep (located above the s band that 
has fallen below Ep) have amplitudes within the a 
cage. We identify these as Px, Py, and Pz bands, re- 
spectively, as confirmed by a fit of the dispersions to 
the tight-binding model, where the fit is again ex- 
cellent. To be precise Ep intersects the Px and Pz 
bands, which are degenerate, reflecting the symme- 
try of the cluster. The fitted hopping integrals for 
s, Px-, Py and Pz bands are, respectively, {tx,ty,tz) — 
(-30.0, -25.0, -0.5), (125, -25.0, -62.5), (-12.5, 25.0, -50.0) 
and (-12.5,-0.5,75.0) in meV. The hopping integrals 
are almost an order of magnitude greater than those for 
KiLTA, as expected from the larger cluster size. An 
ESR experiment shows that the g value decreases for 
K„LTA with n > 2 [13]. This may be understood as 
the degeneracy of the p bands enhancing the spin-orbit 
interaction and reducing the g value. 

Now we come to a big question of whether the sys- 
tem is strongly correlated. As mentioned, a puzzle is 
that experimentally KiLTA and K3LTA are insulators, 
while the LDA finds them to be metals. So we first 
estimate the Coulomb matrix elements. The largest 
one is the intra-orbital Coulomb interaction U, which is 
U — J |(/)i(x)|^y(x - y)\(j)i{y)\^dxdy, where (/>'s are the 
wave functions at T. Here V is assumed to be the bare 
Coulomb interaction as a first approximation, since the s 
and p wave functions are well localized in the interior of 
each cage [14]. U is calculated to be (7 ~ 4.5 eV for the s 
band in KiLTA, and ~ 4.0 eV for the px band in K3LTA. 
Given that U/W ~ 10 > 1, where W is the band width, 
we can expect that these materials arc Mott insulators. 

However, since the relevant bands are p bands with 
very anisotropic dispersions, we have to be careful in es- 
timating the critical Uc for the metal-insulator transi- 
tion. Here we have employed the dynamical mean-field 
theory [15] with the maximum entropy method [16] to 
estimate the transition point by calculating the spectral 
function in the single-band Hubbard model [17] for a typ- 
ically anisotropic {tx : ty : tz =5:1:1) case as well as the 
isotropic (1:1:1) one. 

Fig. 4 shows the spectra for various values of Hubbard 
U/W, where we can see that the system becomes an in- 
sulator (as identified from a gap in the spectral function) 
for U/W larger than ~ 2 [18] in the anisotropic case. 
So we conclude that this particular K-doped zeolite is 
on the Mott-insulator side. In more general terms we 
expect that the metal-insulator transition can be con- 
trollable through control of U/W. For example, a zeo- 
lite called faujasite is known to be metallic when alkali- 
metal doped [19], where this form of zeolite has a signif- 
icantly wider (7A against 5 A for LTA) window between 
the cages, which should result in a smaller U/W. 

We finally come to the magnetic property. To discuss 
this we require the inter-orbital U' and exchange integrals 



J as well, where U' = J |(/)i(x)|2y(x - y)|(/.2(y)|2dxdy, 
J = J 4(x)02(y)^(x - y)0i(y)02(x)dxdy. For the p^^- 
pair in K3LTA, we find U' ~ 3.7 eV, and J ~ 0.7 eV, 
which are similar to those roughly estimated by Nozuc et 
al [4]. Multi-orbital systems are in general favorable for 
ferromagnetism, since the inter-orbital kinetic-exchange 
coupling, Jp"*''"' = —2t^/{U' — J), is ferromagnetic (ac- 
companied by an orbital superlattice structure). Jp"*®'' is 
estimated here to be (40 ~ 100) K for p^ and Pz bands in 
K3LTA. This coupling competes with the intra-orbital 
kinetic-exchange coupling, f^^'^'^ ~ AUt^/{U'^ - J^), 
which is antiferromagnetic. These exchange energies are 
(60 ^ 160) K for K3LTA, which are two orders of magni- 
tude greater than the exchange energy, At"^ /U c± 0{1) K, 
for the (single-band case of) KiLTA, where the difference 
is mainly due to an order of magnitude difference in t. If 
we look at the magnetic phase diagram in the literature 
[20,21], K3LTA is, as indicated in the inset of Fig.4, right 
around the ferromagnetic phase boundary. 

Experimentally, ferromagnetism with finite magnetiza- 
tion has been observed for K„LTA with 2 < n < 7 [3,4], 
where the Curie temperature 10 K for n ~ 3 while the 
spin susceptibility for T > 30 K exhibits a Curie- Weiss 
law, 1/(T — 6), with the Weiss temperature ~ K for 
KiLTA, while (-40 30) K for K3LTA [3,4]. The dop- 
ing dependence as well as the energy scale are consistent 
with the exchange interactions estimated here. 

Finally we comment on the Mott transition, for which 

K„LTA is experimentally [22] insulating even when the 
nominal doping level (averaged n) is fractional. This may 
possibly be related to a coexistence of differently doped 
regions with a domain structure [22,23]. It would be 
interesting to know whether superconductivity as in the 
high-Tc cuprates can appear when we realize doped Mott 
insulators in the present system (by, e.g., degrading such 
a domain structure) . Future work should also include an 
elaboration of the LDF-I-DMFT approach. These will en- 
able us to systematically study electronic and magnetic 
effects in the "supercrystal" , for which some experimen- 
tal and theoretical attempts are under way. 
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FIG. 1. Left: Geometry of the LTA structure. Right: Unit 
cell (i.e., a cage) of the undoped zeolite considered here, with 
dark blue: Si, light blue: Al, dark green: oxygen, orange and 
red: K. 
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FIG. 2. The band structure of the undoped zeolite LTA. 
The top right panel is a blowup of the lowest conduction band, 
where the red line is a tight-binding fit. Bottom panels are 
wave functions (red contours: positive, blue: negative) at F 
in the bands A-F as labeled in the band structure. 
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FIG. 3. Band structure of a K-doped zeolite K3LTA. Col- 
ored curves are a tight- binding fit. Top left inset depicts the 
atomic configuration, where four red K's (the doped K -I- the 
red K in Fig.l) form a square in the a cage. Bottom panels 
are wave functions at F in the bands A-E as labeled in the 
band structure. 
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FIG. 4. Spectral function, piuo), obtained with the dynami- 
cal mean-field theory for the Hubbard model on the cubic lat- 
tice with the isotropic(a) or an anisotropic(b) dispersions for 
various values of U/W . The bottom inset shows the magnetic 
phase diagram for the degenerate Hubbard model against the 
inter-orbital Coulomb (L'^') and the Hund coupling (J) inter- 
actions after Ref. [21], where the parameter region for the 
doped zeolite is indicated by blue hatch with its widths rep- 
resenting the dependence on the atomic configuration of the 
cluster and on the assumption oi U = U' + 2J on which the 
diagram is drawn. W is here defined as the width of the 
gaussian density of states. 



